Ground State of the Quasi-ID BaVSs resolved by Resonant Magnetic X-ray Scattering 
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Resonant-magnetic x-ray scattering (RMXS) near the vanadium L2,3-absorption edges has been 
used to investigate the low temperature magnetic structure of high quality BaVS3 single crystals. 
Below TV = 31 K, the strong resonance revealed a triple- incommensurate magnetic ordering at wave 
vector (0.226 0.226 £) in the hexagonal notation, with £ = 0.033. The simulations of the experimental 
RMXS spectra with a time-dependent density functional theory indicate an antiferromagnetic order 
with the spins polarized along a in the monoclinic structure. 

PACS numbers: 75.25.-j,71.27.+a,75.50.-y,73.63.-b 



BaVS3 and related compounds have attracted consid- 
erable scientific interest in the last decades, mainly due 
to the presence of both itinerant and localized electronic 
states which involves a subtle interplay between charge, 
orbital and spin degrees of freedom |l| . It shows a very 
rich phase diagram, with a metal-insulator (MI) transi- 
tion at Tmi = 69 K and another transition at around 
Tx — 30 K that involves spin degrees of freedom. 
However, as most of the research on BaVS3 has been done 
on powder samples, information about the low temper- 
ature structure is limited. Nuclear magnetic resonance 
(NMR) and nuclear quadrupole resonance measurements 
revealed a non magnetic ground state below Tx and in- 
terpreted the asymmetric field gradient at the V sites 
as a consequence of orbital ordering [2|. Susceptibility 
measurements also indicated a spin pairing ground state 
below Tx [H, while heat capacity measurements did not 
detect any transition Q ■ The development of an orbital 
ordering between Tmi and Tx has also been suggested 
from dielectric function measurements On the other 
hand, powder neutron diffraction revealed very few mag- 
netic Bragg peaks sustaining a magnetic ground state |6j . 
Supposing a ferromagnetic (FM) coupling along c, the 
structure refinement suggested an incommensurate long- 
range ordering with the wave vector (0.226 0.226 0)h 
(hexagonal lattice coordinates) in reciprocal lattice units 
(rlu). A magnetic origin of this transition has also been 
supposed to explain muon spin resonance measurements 
[11]. Therefore the precise nature of the low temperature 
structure needs to be clarified. 

In order to determine the exact nature of the low tem- 
perature structure, we used soft resonant magnetic x-ray 
scattering (RMXS), which has proved to be an excellent 
tool for probing magnetic ordering in 3c? materials 0- 



Hoj . The large resonant enhancement of the scattering 
cross section, in combination with the high brilliance of 
the x-ray beam provided by third generation synchrotron 
facilities enable the investigation of small single crystals. 

In this Letter we report the direct observation of an 
antiferromagnetic (AFM) ground state, by RMXS at the 
V-L edge. Below Tn — 31 K, the strong resonance 
clearly shows an in-plane modulation, which agrees with 
a previous suggestion, and brings to light a small incom- 
mensurate component in the c-direction. Through time- 
dependent density functional theory (TD-DFT) calcula- 
tions, we reproduce the experimental data and show that 
the spins are polarized along the a direction in the mon- 
oclinic structure. 

At room temperature the S = 1/2 BaVS3 (3d 1 ) sys- 
tem is metallic and its PQ^/mmc structure consists of 
a hexagonal packing of quasi ID-chains of face-sharing 
VS6 octahedra along c, with two V atoms per unit cell 



(uc) A zigzag deformation of the V-S chains at 



240 K doubles the number of (in-phase) chains in the uc 
and reduces the crystal symmetry to an orthorhombic 
structure (Cmc2i) (H-0 transition). This transition is 
accompanied by the creation of 6 twin domains in the 
aft-plane [T^. Peierls instability drives the system to an 
insulating phase at Tmj, lowering further the symmetry 
to the Im space group [131 ]. The 2k p charge density wave 
(CDW) is accompanied by a tetramerization of the V 
chains. At Tmi the magnetic susceptibility shows an an- 
tiferromagnetic (AFM)-like drop, indicating that the spin 
degrees of freedom are strongly affected by the MI tran- 
sition [l,[l3]. On the other hand, the critical wave vector 
2kp = 0.5c* (c* stands for the reciprocal lattice vector 
related to c) indicates that only one of the two d electrons 
per unit cell participates directly in the CDW [3]. The 
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FIG. 1: (Color online) Energy dependence of the magnetic 
reflection at (0.226 0.226 £) measured at 10 K. The solid 
line is the result of the FDMNES code calculation (see text). 
The inset shows the V-I/2,3 x-ray absorption spectra recorded 
in total fluorescence yield mode (dots) and compared with 
FDMNES simulation (solid line). 



effective moment deduced from the Curie- Weiss behavior 
of the magnetic susceptibility is 1.223 hb, consistent with 
one localized spin every other V site [13J . Band structure 
calculations reveal two di fferent t 2g active electron states 

two narrow E g bands and one 



at the Fermi level 
dispersive A\ g band extending along c undergoing the 
CDW instability [13]. 

The experiments were performed at the UE46-PGM1 
beamline at the BESSY II synchrotron light source in 
Berlin, using a two-circle diffractometer with horizon- 
tal scattering geometry. The sample was mounted on a 
copper goniometer and attached to a continuous He-flow 
cryostat cooling the sample down to 10 K. The signal 
was recorded with a photodiode with a vertical (hori- 
zontal) acceptance of 4° (0.4°) providing a momentum- 
resolution of 0.001 rlu. The BaVS3 single crystal, grown 
as explained elsewhere [l8j], with a size of 2 x 3 x 0.5 mm 3 , 
was aligned along [110]//, allowing a rotation about the 
scattering vector. The results shown here are mainly ob- 
tained on an uncleaved surface and subsequently con- 
firmed on an in-situ cleaved surface. 

At 10 K, a strong resonant enhancement of 
the scattered intensity is observed at wave vector 
(0.226±0.001 0.226±0.001 £) H near the V-L 2:3 absorp- 
tion edges (Fig. [T]). This result, is the first direct ob- 
servation on a single crystal, of the in-plane components 
of the magnetic ordering vector suggested from powder 
neutron diffraction [f| . The x-ray absorption spectra (in- 
set Fig. [1]) shows two main structures with maxima at 
517 eV and 524 eV, corresponding respectively to the 
V-L3 (2p 3 / 2 — > 3d) and V-L2 (2pi/2 —> 3d) dipole tran- 
sitions. The resonant signal presents two double peak 
structures near 515 eV and 523 eV close to the V-L3 



and L 2 absorption edge, respectively (Fig. [T]). It can be 
attributed to different t 2g states and related to two tem- 
perature dependent structures observed by total electron 
yield measurement (l9| . This interpretation is also sus- 
tained by the calculated density of states (not shown) 
which shows an admixture of several t 2g states at the 
Fermi level. Thereby the RMXS signal probes mainly 
the magnetic moment of the t 2g states while the e g states 
make a negligible contribution. The scattered intensity 
disappears at Tn — 31±1 K (see inset Fig. [2]) at the 
same transition temperature reported in previous exper- 
iments d, @, |2(j. Moreover, the intensity related to the 
square of the order parameter decreases nearly linearly 
when the temperature increases, as expected from a sec- 
ond order transition within the mean field approxima- 
tion. 

To determine precisely £, the rocking curve of the mag- 
netic signal was measured at two different azimuthal an- 
gles * = 0° and * = 90° (Fig.[2|). The scattering geome- 
try and the definition of are illustrated in Fig. [3l Note 
that at \& = 0° (90°) a and b (c) are in the scattering 
plane. At ^ = 0° several structures related to differ- 
ent twin domains are clearly observed, while at \t = 90° 
the rocking curve shows only two well separated peaks. 
These two peaks, separated by A9 = 4.77°, are a fin- 
gerprint of two incommensurate magnetic reflections at 
(0.266 0.226 ±£). When * = 90° the magnetic signal 
from different twin domains is integrated vertically and 
can not explain the two structures. A9 is therefore a 



y = 90° 




54 56 
e (deg.) 

FIG. 2: (Color online) Rocking curve of the magnetic signal 
measured at ^ = 90° (upper panel) and at $ = 0° (lower 
panel). The blue lines correspond to a fit with a Gaussian 
profile. The inset shows the temperature dependence of the 
(0.226 0.226 0.033) magnetic reflection measured at 10 K with 
516.5 eV photon energy. 
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FIG. 3: (Color online) Azimuthal-angle (ty) dependence of the 
ratio of resonant-peak intensities measured with a and n in- 
cident polarized beam (black triangles). The blue dashed line 
corresponds to the FDMNES simulation when the spins are 
polarized along a/ m (see text). The schematic illustrates the 
definition of $ and shows the scattering geometry at $ = 0°. 



parameter able to determine accurately the incommen- 
surability along c*, with a calculated value £ = 0.033 
± 0.001 in rlu 
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Its small value indicates that the 
interactions are mainly FM in agreement with previous 
predictions. 

In order to determine the spin polarization, the mag- 
netic signal over a complete azimuthal-period has been 
recorded on one twin-domain. To avoid surface and sam- 
ple size effects, Fig.[3]shows the ratio Ia/I^ where I a and 
I v are the intensities recorded with incoming polarization 
respectively perpendicular and parallel to the scattering 
plane. The spectra have a period 7r with a maximum 
at * = 180° and a minimum at * = 90°. Due to the 
complexity of the structure and the multi-electronic phe- 
nomena occurring in the resonant scattering process, we 
used simulations to determine the best magnetic struc- 
ture compatible with the experiments. 

Simulations were performed with the FDMNES code 
[22j |. It has already been extensively used to simulate 
resonant elastic scattering spectra in a large variety of 
compounds, mainly at isT-edges. To take into account the 
strong photoelectron-hole interaction occurring at the 
-^2,3-edges of transition metals, time-dependent density 
functional theory (TD-DFT), was recently incorporated 
in the code [23[ . It was first developed by Runge et al. [24 1 
and later applied to simulate x-ray absorption spectra in 
3d transition metal compounds [25] . We present here its 
first use in the context of RMXS. Calculations are fully 
relativistic and use a local exchange correlation Kernel 
(so called TD-LSDA) including the spin-orbit coupling 
for both the core and valence states without any ap- 



proximation. To calculate the resonant atomic scattering 
amplitude, we use the multiple scattering frame in 5 A 
radius clusters around each absorbing atom containing 
around 23 atoms. 

Starting from the low temperature I m monoclinic cell 
based on the crystal structure of Fagot et al. [13] , several 
magnetic models have been considered. The simulations 
were done in the simplest 1x2x1 monoclinic supercell, 
with the spins polarized in the ab plane and an AFM 
coupling along bi m (Fig. |4}. In this supercell 4 layers of 

V atoms contains 8 inequivalent V. One layer is schemat- 
ically represented in Fig. 01 where the 2 non-magnetic V 
atoms are located at the nodes of the spin wave. No spin 
frustration takes place and one spin S = 1/2 every two 

V atoms in the uc is in agreement with Refs. [ID, [r| [26[ . 
The calculations are done with the (0,-1/2, 0)/ m mag- 
netic reflection equivalent to (h' ,-2h' ,0) h , with h' = 0.25. 
Due to the 6-fold symmetry of the hexagonal structure 
this reflexion is experimentally indistinguishable from the 
(h\h',0)H- The calculated XAS and resonant scattering 
spectra reproduce the overall shape of the experimental 
data (Fig. [T]) and the excellent agreement of the peak 
positions indicates the correct treatment of the multi- 
electronic effects. 

Simulations on a larger supercell, like the 1x20x1 
with a 206/ m /9 modulation (/i' = 0.225 very close to 
h = 0.226, the experimental indexation), does not 
change the results of the calculation. Adding a mag- 
netic moment on the central V atoms does not influence 
the shape of the spectra but introduces a modification 
in the intensity. The calculation in the 4x4x1 supercell, 
corresponding to the (h',h\0)H reflection, leads to a spin 
wave along (110)# and does not affect the simulated az- 
imuthal dependence. Moreover this cell is less reliable 
because the AFM ordering extends along the ai m +bj m 
direction where the V-V distances are longer. 

The best I a /I-K simulation is obtained when the spins 
on the V atoms are polarized along aj m (Fig. [3]). A 
clear maximum (minimum) is seen at ^ — 180° (\& = 
90°) in agreement with the experiment. This result is 
also supported by the susceptibility (\) measurements 
which pointed out that below TN,_Yab is smaller than \c, 
indicating an easy axis within ab [3j . When the spins are 
polarized along bj m the simulated azimuthal spectra is 
shifted by a 7r/2 phase with respect to the experiment, 
and with a magnetic moment oriented along c/ m , I a / I-n 
is equal to unity. 

Usually the RMXS data are analyzed with the stan- 
dard formula |(e* x e s ) • rh\ 2 [13,[2i|, where q (e s ) is the 
polarization of the incoming (scattered) photons and rh 
the unit vector of the density of magnetic moment. In 
spite of the importance of the local symmetry in the de- 
termination of magnetic structures, not taken into ac- 
count in this formula [53, [5^] , it is valid here due to the 
orientation along a high symmetry direction of the mag- 
netic moment. In fact, the simulation of the azimuthal- 
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FIG. 4: (Color online) Magnetic structure of BaVS3 in the 
monoclinic (black) and hexagonal (blue) lattice deduced from 
combined RMXS measurements and ab inito calculations. 

angle dependence with the above mentioned formula also 
indicates that the spins are oriented along aj m . This 
attests that the multi-electronic effects do not strongly 
affect the sensitivity to the direction of the magnetic mo- 
ment and shows that, the density of magnetic moment 
remains parallel to the total magnetic atomic moment. 

In summary the strong resonant enhancement ob- 
served below Tjv = 31 K at the V-L2.3 edge at reflec- 
tion (0.226 0.226 0.033)//, reveals a triple magnetic in- 
commensurate ordering in BaVS3. The temperature and 
azimuthal-angle dependence of the signal attests to the 
magnetic origin of the transition at Tjy. Through simu- 
lation of the resonance with the FDMNES code in a TD- 
DFT multiple scattering frame, we reproduced the exper- 
imental data and established that the spins are polarized 
along aim- The long-range modulation along c shows 
that the interactions are predominantly FM, indicating 
either a conical structure with the magnetic moment 
slightly deviated from the a-axis or a spin-density wave 
along c. Further theoretical investigations are needed to 
explain the precise origin of the out-of-plane modulation. 
We believe that the determination of the magnetic struc- 
ture in BaVS3 will help ongoing efforts for further theo- 
retical work which may help to understand the puzzling 
interplay between magnetic, charge and orbital order in 
low-dimensional systems. 
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